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I 


INTRODUCTION 


Considerable  Interest  has  been  generated  in  using  Fourier  transform 
matched  filters  for  object  discrimination  and  tracking.  In  order  to  recognize 
several  objects  or  several  orientations  of  the  same  object  many  matched 
filters  may  be  needed.  Interest  has  been  shown  in  storing  several  of  these 
filters  in  a  small  area  [1,2].  This  has  been  done  by  multiply  exposing  the 
photographic  plate  or  by  using  arrays  produced  by  holographic  elements  or 
other  means  [3,4].  Arguments  are  presented  here  which  investigate  how  closely 
packed  these  filters  can  be  and  still  have  good  resolution  between  them.  In 
this  instance  close  packing  refers  to  double-exposing  a  photographic  plate  to 
an  input  scene  at  two  slightly  different  rotation  positions.  Classic  line- 
shapes  are  used  to  model  examples  of  experimental  results  and  some  conclusions 
are  drawn  using  the  Rayleigh  criteria  as  a  guideline. 


II 


THEORY 


The  Gaussian  intensity  distribution  has  the  form 
li  “  I0  e"2r2/w2 


(1) 


where  r  is  the  radius  from  the  center  of  the  distribution  and  2w  is  the  full 
line  width  at  I  =(l/e2)l0-  If  another  distribution  of  the  same  form  is 
located  a  small  distance  away 

I2  -  I0  e-2<r-6l)2/w2  (2) 


then  the  sum  of  these  distributions  is 


Let 


lo  «’ 


2r2/w2 


-2(r-6j)2/ 


w^ 


I  (r=»0) 

I (r=  fij/2) 


Substituting  from  the  above  gives 


P 


1 


1/2 


+ 


(3) 


(4) 


(5) 


PI  is  thus  the  ratio  of  the  peak  at  r  *  0  to  the  height  of  the  valley  between 
the  two  distributions.  Solving  for  using  the  fact  that  the  second  term  is 
much  smaller  than  the  first  gives 
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w 


2  Ln  (2p) 


(6) 


Using  the  Rayleigh  criteria  as  a  rough  guideline,  pj  ■  1.23  and  substituting 
this  into  Eq.  (6)  yields 


*  1.34  w 


(7) 


The  same  sort  of  analysis  may  be  applied  to  a  Lorentzian  distribution  having 
the  form 

j  =  _ A _  +  _ A _ 

r2  +  (o/2) 2  (r-62)2  +  (a/2)2  (8) 

where  A  is  a  constant  equal  to  I0(c/2)2,  r  is  the  radius  and  a  the  full 
linewidth  at  half  the  maximum  intensity.  This  produces  a  62  of 


62  =  o/2  j^4p2  -  3  16p2(p2-i)+i  j  ^2 

Then  for  P2  a  1.23  as  a  value  of  the  intensity  ratio 
&2  *  1.03  a 

For  a  sine2  distribution 


/  \  2 

sin  a  (r-<5_) 

Is  in  a  r  1  T 

3 

V  ”  )  ° 

a/r-5^) 

(9) 


(10) 
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where  again  r  is  the  radius  from  the  center  of  the  distribution  and 

a  =  u/r(o)  (12) 

where  r(o)  is  the  value  of  r  at  the  fist  zero  of  intensity.  This  distribution 
yields 


p3 


a2  6^ 2  +  sin2  6^ 
8  sin2  (a6^/2) 


(13) 


which  must  be  solved  graphically  or  by  iteration  using  p3  =  1.23.  For  small 
values  of  63  ,  note  that  a2  63  2  >  >  sin2^. 


III.  EXPERIMENTAL  RESULTS 


Figure  1  shows  three  distributions  compared  with  experimental  data  taken 
using  visible  (Helium-Neon)  real-time  optical  correlation  methods.  The  input 
scene  is  fed  via  an  RCA  television  camera  (model  number  TC  1005)  into  a 
Videotek  Monitor  (model  VM-12PR)  which  is  used  as  the  input  to  a  Hughes 
liquid  crystal  light  valve  (the  same  one  used  in  Ref.  4)  which  produces  the 
coherent  image  that  is  fed  into  a  real-time  Vander  Lugt  type  correlator. 
Experiments  of  this  type  are  described  in  detail  elsewhere  [5].  All  three 
distributions  have  been  arbitrarily  fit  to  experimental  data  at  r  =  0  and  r  =* 
2.94°  using  w  =  2.80°,  a  *  2.07°,  and  a  =  .775  deg~l.  This  gives  6j  =  3.8°, 

62  **  2.13°,  and  63  *  4.0°.  Note  that  these  widths  should  depend  upon  the 

spatial  frequency  of  the  input  scene  as  well  as  other  factors  and  thus  serve 
only  as  examples  here.  In  this  particular  experiment  the  input  scene  was  a 
low  resolution  black  and  white  aerial  photograph  of  Huntsville,  Alabama.  The 
photograph  was  placed  on  a  rotatable  stage  and  matched  filters  made  at  the 
desired  rotation  angles.  Experimental  data  in  Figure  1  is  the  correlation 
intensity  for  a  filter  made  at  0°  rotation  of  the  input  scene.  After  develop¬ 
ment  of  the  film  plate  (Kodak  649F)  it  was  replaced  in  the  correlator  and  the 
input  scene  rotated  slowly  as  the  corresponding  correlation  intensity  was 
measured  using  an  RCA  CCD  camera  fed  into  a  television  monitor  then  digitized 
by  a  Colorado  Video  Analyzer  (model  321).  This  signal  was  then  recorded  by  a 
Hewlett-Packard  model  680  strip  chart  recorder.  Figure  2  is  representative  of 
actual  data  taken  as  the  input  scene  was  rotated  through  the  position  (0° 
rotation)  where  the  matched  filter  was  made. 


In  order  to  investigate  multiply  stored  filters,  a  single  plate  was 
exposed  twice  with  different  rotation  angles  of  the  input  scene.  After  deve¬ 
lopment  using  standard  techniques  for  Kodak  649F  plates,  the  filter  was  rein¬ 
serted  into  the  correlator  and  data  taken  as  before.  Figure  3  shows  the 
actual  data  and  a  plot  of  Equation  8  using  the  Rayleigh  criteria.  The  expo¬ 
sures  were  made  using  the  same  criteria  for  62  of  2.13°  of  scene  rotation. 

The  agreement  is  reasonably  good.  Some  variation  is  expected  due  to  the  dif¬ 
ficulty  in  obtaining  equally  intense  correlations  for  the  two  filters.  This 
depends  upon  the  exposure  times  used,  film  response,  and  several  other 
variables.  In  the  theoretical  analysis  the  calculations  could  easily  have 
been  done  using  distributions  having  different  peak  heights  but  then  the  value 
of  6  would  depend  upon  these  heights  which  can’t  be  known  with  any  certainty 
in  advance  of  doing  the  experiment.  In  the  calculations  presented,  it  was 
also  assumed  that  individual  correlation  intensities  should  be  added  directly 
rather  than  adding  the  corresponding  fields  then  squaring  for  the  intensity. 
This  has  been  justified  experimentally  by  storing  one  filter  made  of  0°  scene 
rotation  and  another  filter  at  90°  scene  rotation  then  adding  the  resulting 
correlation  intensities  algebraically.  This  distribution  was  then  compared 
with  data  similar  to  Figure  3.  The  results  were  found  to  be  essentially 
equal.  This  indicates  that  there  is  little  if  any  phase  addition  contribution 
to  the  sum  of  the  two  correlation  intensities. 


IV.  CONCLUSIONS 


In  this  brief  report  a  method  of  analyzing  closely  spaced  optical  matched 
filters  has  been  presented  along  with  experimental  data  which  tends  to  support 
the  findings.  The  use  of  the  Rayleigh  criteria  may  not  be  strictly  correct  in 
that  this  criteria  is  for  a  Bessel  function  intensity  distribution  only  but  it 
serves  as  an  initial  guideline  [6].  Experimental  data  seems  to  suggest  that  the 
separation  criteria  will  become  more  stringent  as  laboratory  methods  improve. 


SCENE  ROTATION  (DEGREES) 


Intensity  (arbitrary  units) 


Correlation  Intensity  as  Input  scene  was  rotated  through  0°  where  matched  filter  was  made 
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